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After reviewing some basic concepts of the theory of strongly interacting 
matter above normal nuclear energy density and reviewing some salient 
results of the experimental program at the Rclativistic Heavy Ion Collider 
(RHIC), these lectures explain why the quark-gluon plasma observed in 
the RHIC experiments has been called a "perfect liquid." They then give 
an introduction to some recent ideas concerning the possible origin of the 
nearly inviscid nature of the quark-gluon plasma and discuss the connection 
between low viscosity and strong parton energy loss of hot QCD matter. 

PACS numbers: 

1. QCD matter under extreme conditions 

Nuclear matter can be compressed in two very different ways: A slow 
squeeze results in cold matter at high net baryon density; a rapid squeeze 
produces hot matter high in energy density. The first type of compression is 
impossible to perform in the laboratory, but it can be achieved in the cosmic 
environment by stellar collapse: The cores of old neutron stars are thought 
to contain nuclear matter at up to ten times normal baryon density and at 
rather low temperatures on the nuclear scale, in the keV range. Hot, dense 
nuclear matter can be produced in the laboratory by colliding two heavy 
nuclei at high energies. The higher the collision energy, the lower is the net 
baryon content of the matter formed in the center of mass. Of course, hot 
nuclear matter also existed for a brief time shortly after the Big Bang, with 
temperatures exceeding 200 MeV (about 2 x 10 12 K) for the first 10 /js in 
the history of our universe. 

Nuclear collisions of this kind occur regularly in cosmic ray reactions, 
but these are rare and difficult to study. Systematic studies using high- 
energy accelerators now have a history of over 30 years, beginning with 
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the BEVALAC at Berkeley, followed by the BNL-AGS and the CERN-SPS, 
and culminating presently at the Relativistic Heavy Ion Collider (RHIC) at 
BNL. In the near future, even higher energies will be explored at the CERN- 
LHC, and large net baryon densities will be studied at RHIC and, later, at 
GSI-FAIR. The goal of all these experimental programs is to produce and 
study equilibrated nuclear matter at energy densities far exceeding that of 
ground state nuclear matter (so ~ 0.14 GeV/fm 3 ). As we shall discuss 
below, this goal has been achieved with spectacular results, especially at 
RHIC. 

The theory of nuclear matter under extreme conditions is based on 
the fundamental theory of strongly interacting matter, quantum chromody- 
namics (QCD). This gauge theory of SU(3)-color describes the interactions 
among quarks and gluon by means of the Lagrangian 

Lqcd = -\G%G a ^ + ]T (id, - gAp a ) ^ + J2 m f^f- (!) 

/ / 

In this expression, A® denote the gluon potential and field strength 
tensor, respectively, with the color index a running from 1 to N% — 1 = 8. 
ifjf denotes the quark field of flavor / = u,d, s, . . ., 7^ are the Dirac matrices, 
t a the generators of SU(3) in the fundamental representation, and nif the 
(current) quark masses. 

At high density, the separation between quanta becomes small and either 
the Fermi momentum or the thermal momentum becomes large. The asymp- 
totic freedom of QCD then suggests that interactions among the quarks and 
gluons are relatively weak, and their contribution to the energy density of 
the matter should be relatively small. We can then estimate the energy den- 
sity as a function of the temperature T. Denoting the number of internal 
degrees of freedom by v and assuming zero net baryon density (hb = 0) we 
have 

where we have neglected the quark masses in the last step. The degrees of 
gluons and quarks are easily counted: vq = 2 (A,? — 1) [spin and color] and 
vq = 4N c Nf [spin, baryon number, color, and flavor]. For Nf = 2 flavors 
(u, d) one finds that a temperature of 160 MeV is needed to reach an energy 
density of 1 GeV/fm 3 , roughly 7.5 times that of normal nuclear matter. 

If one wants to obtain more reliable information about the equation of 
state of QCD matter or wants to find out in which temperature domain 
the rough estimate ([2]) applies, one must turn to exact calculations of the 
energy density. At present, this requires massive computer simulations of 
QCD discretized on a lattice. Over the past few years, increasingly precise 
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lattice calculations of thermal QCD, extrapolated to the continuum and 
thermodynamic limits and to small quark masses, have become available 
(see Fig. [TJ. These show a rapid transition from a low-energy region to a 
region with nearly constant e/T A at T c ~ 160 MeV. Detailed studies of the 
baryon number susceptibilities have shown that below T c the matter can be 
understood as a rather dilute gas of hadrons. Above T c , the matter behaves 
as a SU(3)-colored plasma of interacting quarks and gluons. The various 
susceptibilities (specific heat, scalar quark density, Polyakov loop) exhibit 
strong peaks at T c , but they do not diverge. The transition between the 
low-temperature and the high-temperature domains is thus not a true phase 
transition, but has the character of a rapid cross-over pQ . Precisely how wide 
the cross-over region is and whether it is characterized by a universal value 
of T c or slightly different values for different physical observables, is a matter 
of intense discussion j2j 13] . 
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Fig. 1. Thermodynamic quantities e, P,s divided by T for QCD matter as function 
of temperature T for fiB — 0. The quantity (45/27T 2 )s/T A can be considered as a 
measure of the number of thcrmodynamically active degrees of freedom. 



2. The quark-gluon plasma 



What is clear from the lattice results is that above T c many quantities are 
well described by positing that the thermodynamically active constituents of 
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QCD matter have the quantum number of quarks, not of hadrons. A striking 
example are the off-diagonal flavor susceptibilities. If a quasi-particle picture 
applies to the quark-gluon plasma, the ratio g] 

_ (XY)-(X)(Y) 

of the covariance of the flavor quantum numbers X and Y and the variance 
of Y alone is a sensitive indicator of the flavor quantum numbers of the 
plasma constituents [5]. In the narrow temperature range T = (1 ± 0.1)T C 
lattice calculations show a striking transition of the quantities Cbs an d Cqs 
from the values expected for a hadron gas containing strange baryons with 
integer electric charge Q and baryon number B to matter whose carriers of 
baryon number and electric charge have the fractional quantum numbers 
characteristic of strange quarks (see Fig. [2j . 
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Fig. 2. Flavor susceptibility ratios Cbs an d Cqs as a function of T/T c . The results 
of lattice simulations [B] are compared with expectations for a weakly coupled 
hadronic gas and a quasi-particulate quark-gluon plasma [7]. The transition to 
Cbs — Cqs = 1 expected for matter in which baryon number is carried by quarks 
occurs in a narrow range around T c . 



How high in temperature one needs to go for gluons to also be described 
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in such a picture is still unclear, but the thermodynamic functions are repro- 
duced well in an approach based on perturbative quasi-particles for T > 3T C 
[8j [9] . These gluonic quasi-particle modes acquire a thermal "mass" of order 
gT, which helps avoid many infrared problems of the bare thermal gauge 
theory. An immediate consequence is the screening of the long-range color 
force responsible for quark confinement at zero temperature. In the weak 
coupling limit, the static color-electric potential Aq of a heavy color charge 
embedded in the quark-gluon plasma satisfies the Poisson equation 

-V 2 A a = gp a G (A a ) + gp a Q (A a ), (4) 

where Pq^q(Aq) = —/j,q/qAq are the induced color polarization densities of 
gluons and quarks in the plasma, respectively. The perturbative calculation 
yields 

p 2 G = (gT) 2 ; p 2 Q = N f (gT) 2 /6, (5) 

implying an exponentially screened static color potential with a range Mg/q < - 
0.5 fm. This picture suggested by a perturbative treatment is borne out by 
lattice simulations of the free energy of a pair of heavy, static quarks |10] 
(see Fig. [|). 

For a long time it was thought that extended, though screened, color 
fields could not be generated in a quark-gluon plasma, because there exists 
no coherent source for such fields. (Owing to the absence of free quarks 
in nature, one cannot construct a "color battery.") Over the past decade, 
however, it has been realized that dynamic plasma mechanisms permit the 
formation of extended color fields by a kind of instability driven "color 
dynamo." IT2] These instabilities are generalizations of the phenomena 
first identified by Weibel in electromagnetic plasmas with an anisotropic 
momentum distribution [13]. Fluctuations in the (color) current densities 
of such a plasma are found to induce (color-)magnetic fields which, in turn, 
amplify the original fluctuations producing a run-away effect. Momentum 
anisotropics are expected in a heavy-ion collision, where the longitudinal 
motion along the beam axis provides for a natural source of momentum 
space anisotropy. In a quark-gluon plasma, the exponentially growing fields 
eventually saturate due to the nonlinear interactions present in the Yang- 
Mills equation governing the SU(3) gauge field p3]. Detailed simulations 
(see Fig. |4j) have shown that the fields in the steady state are characterized 
by an intensity of order (E 2 ) ~ (B 2 ) ~ g 2 T 4 and a correlation length of 
order (gT) -1 p2]. 

A second conceptually important difference between normal nuclear or 
hadronic matter and the quark-gluon plasma phase is the restoration of 
the spontaneously broken chiral symmetry due to the disappearance of the 
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Fig. 3. Free energy F\ of a pair of heavy quarks forming a color singlet as a function 
of QQ separation r for various different temperatures. The potential is increasingly 
screened as the temperature increases. The screening effect at T < T c is due to 
pair production of light quarks from the vacuum when R is large. 



vacuum quark condensate (0|gg|0) at high temperature. As a result, the dy- 
namical ("constituent") masses of the light quarks u, d (s) of order 300 MeV 
(500 MeV) give way to the much smaller current masses of order 5 — 10 MeV 
(100 MeV) believed to be generated by the interaction with the Higgs field, 
as illustrated in Fig. [5] It should be noted, however, that the interaction of 
the quarks with the thermal medium induces a new type of dynamical mass 
of order gT, which does not violate chiral symmetry. Although the phe- 
nomenon of chiral symmetry restoration is well established owing to lattice 
simulations, it is not easy to identify experimentally accessible signatures 
of this phenomenon. Electromagnetic probes, such as lepton pairs, may 
provide the best handle, but it remains uncertain how the effects of quark 
deconfinement and chiral symmetry restoration in the emission spectrum 
can be separated. 

Rigorous algorithms for the simulation of lattice QCD only exist for zero 
net baryon density. In recent years, however, several ingenious techniques 
have been invented that allow to extrapolate the simulations to nonzero 
values of the baryon chemical potential \x b [HI US Qj2 EQ] . Most of these 
simulations predict that the cross-over regime between hadronic matter and 
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Fig. 4. Spatial and temporal structure of dynamically generated color fields in 
a spatially one-dimensional quark-gluon plasma with an anisotropic momentum 
distribution |16j . 




a s c 6 ' condensate 

Fig. 5. Right panel: The overwhelming part of the mass of "constituent" u, d, and 
s quarks in hadrons in the normal QCD vacuum are dynamically generated by 
the vacuum quark condensate (bottom). This contribution to the mass disappears 
above T c and only the current mass generated by the Higgs field (top) remains. 
Left panel: Current mass due to Higgs field (blue) and constituent mass due to 
quark condensate (red) of the six known quark flavors. The quark mass is shown 
in MeV on a logarithmic scale. 
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quark-gluon plasma narrows as ^b increases and eventually gives way to 
a discontinuous phase transition. The critical endpoint of the first-order 
transition line may lie as low as fiB ~ 300 MeV [21 . Unfortunately, the 
precise location of the critical point has been found to be extremely sensitive 
to the values of the current quark masses. Some calculations even indicate 
that the critical point might not exist [22]. Following the phase boundary 
to even higher values of \xb one may encounter a triple point, where the line 
separating hadronic matter from the quark-gluon plasma meets the phase 
boundary between quark-gluon plasma and color superconducting quark 
matter |23j . It must be noted, however, that such studies are still very 
much exploratory and speculative. 

3. Results from RHIC 

The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory was contructed to heat nuclear matter beyond T c by colliding 
heavy nuclei at sufficiently high energy, and thus to discover the quark-gluon 
plasma and explore its properties [23]. Because lattice QCD simulations 
are not yet able to reliably address the dynamic processes that govern the 
plasma evolution and response to perturbations, Initial expectations for the 
dynamic properties of the quark-gluon plasma were based on approximate 
analytical QCD calculations which assume that the color forces between 
quarks and gluons at high temperature are weak and can be treated per- 
turbatively. Such analytical calculations suggested that the quark-gluon 
plasma should behave like a dilute gas, a loose assembly of particles which 
explodes in a more or less spherical pattern. The RHIC experiments have 
provided spectacular evidence that the quark-gluon plasma defies this ex- 
pectation, at least in the temperature domain T < 2T C accessible at RHIC. 
Far from behaving like a dilute gas, the matter created in collisions of heavy 
nuclei flows like a nearly "perfect" liquid of liberated quarks, with almost 
undetectable viscosity. 

The insights into the nature of quark-gluon plasma have made the first 
phase of RHIC a great success. In its first six runs (2000-2006), RHIC 
has provided four different collisions systems (Au+Au, d+Au, Cu+Cu, and 
p+p) at a variety of energies, ranging from nucleon-nucleon center-of-mass 
energies of 19.6 to 200 GeV. The largest data samples were collected at the 
highest energy of 200 GeV, where the accelerator has achieved sustained 
operation at four times the design luminosity. The ability to study proton- 
proton, deuteron-nucleus, and nucleus-nucleus collisions at identical center- 
of-mass energies with the same detectors has been key to systematic control 
of the measurements. The experiments also measure the impact parameter 
(distance of closest approach) and the orientation of the reaction plane in 
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each event. This detailed categorization of the collision geometry provides 
a wealth of differential observables, which have proven to be essential for 
precise, quantitative study of the plasma. The results obtained by the four 
RHIC experiments (BRAHMS, PHENIX, PHOBOS, and STAR), published 
in nearly 200 journal articles and four summary publications [25 1 126 | 127 1 128] . 
are in remarkable quantitative agreement with each other. 

The initial set of heavy ion results from RHIC has provided evidence 
for the creation of a new state of thermalized matter at an unprecedented 
energy density of (30 — 100)eo> which exhibits almost ideal hydrodynamic 
behavior. Important results from the RHIC experiments include [29J: 

• chemical (flavor) and thermal equilibration of all observed hadrons 
including multi-strange baryons; only the reaction gg — * ss is known 
to achieve this on the time-scale of the nuclear reaction; 

• strong elliptic flow, indicating early thermalization (at times less than 
1 fm/c) and a very low viscosity of the produced medium; 

• collective flow patterns related to independently flowing valence quarks, 
not hadrons; 

• strong jet quenching, implying a very large parton energy loss in the 
medium and a high color opacity of the produced matter; 

• strong suppression of open heavy flavor mesons at high transverse 
momentum, implying a large energy loss of heavy (c and b) quarks in 
the medium; 

• direct photon emission at high transverse momentum that remains 
unaffected by the medium; 

• charmonium suppression effects that are similar to those observed at 
the lower energies of the CERN-SPS. 

Among the results, three discoveries are most relevant to this lecture: 

"Elliptic" flow: The measured hadron spectra and their angular dis- 
tributions in non-central heavy-ion collisons reveal the enormous collective 
motion of the medium. In addition, measurements of electrons from the 
decays of hadrons containing charm quarks indicate that even heavy quarks 
flow with the bulk medium. These conclusions are based on measurements 
of the second Fourier coefficient v 2 of the azimuthal distribution of particles 
around the beam axis |30[ [31] : 
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where the angle <p p is measured with respect to the reaction plane. The 
RHIC data are in good agreement with predictions for the hydrodynamic 
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expansion of a nearly viscosity-free liquid - often referred to as a "per- 
fect" liquid - from an initial oval-shaped configuration characterized by an 
anisotropic pressure gradient. The relative abundances of the produced 
hadrons with transverse momenta below 1.5 GeV/c, the shapes of their 
transverse momentum spectra, and their anisotropic or elliptic flow patterns 
are well described by ideal relativistic hydrodynamics with an equation of 
state similar to the one predicted by lattice QCD. Ideal hydrodynamics pre- 
dicts a systematic fine structure in the behavior of the elliptic flow strength 
parameter V2 as a function of transverse momentum px for hadrons of differ- 
ent mass (see Fig. [6]). The magnitude of the observed collective flow points 
to rapid thermalization and equilibration of the matter on a time-scale of 
less than 1 fm/c [32 . 
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Fig. 6. The momentum dependence of the flow anisotropy v 2 ("elliptic flow") of 
hadrons measured in on-central Au+Au collisions at RHIC shows the characteristic 
dependence on hadron mass predicted by ideal relativistic hydrodynamics. 

In reality, hydrodynamics is never ideal; any (non-super-)fluid has non- 
vanishing viscosities rj, £, which describe the deviation of the stress-energy 
tensor 

= euVu - PA^ + rj |v/ + W - -tf^V • it) + CA^V • u (7) 

from its ideal form in the presence of flow gradients. The expression ([7| 
is valid in the so-called energy frame and makes use of the abbreviations 
A^ v = g^ u — u^u u and V M = A^ v d v . rj is called the shear viscosity, £ the 
bulk viscosity. £ vanishes for a scale invariant fluid and is thus expected to 
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be small for QCD matter, except possibly in the immediate vicinity of T c , 
because scale invariance is only broken in QCD by the small current quark 
masses and by quantum effects. One of the exciting theoretical discoveries 
of the past few years is the insight that there may be a lower bound on 
the ratio between the shear viscosity r/ and entropy density s of any fluid: 
Airq/s > 1 |33j . A "perfect" liquid is a fluid that attains this lower bound. 
There is mounting evidence from analysis of RHIC data that the matter 
produced is nearly such a perfect liquid, with a viscosity to entropy density 
ratio not larger than a factor of four times the lower bound. 

Quark recombination: Evidence that the medium is composed of de- 
confined, thermalized and collectively flowing quarks comes from detailed 
measurements of the spectra of a wide variety of hadrons. Baryons, which 
contain three valence quarks, show yields that are strongly enhanced rel- 
ative to those of mesons, containing a valence quark and anti-quark, at 
intermediate transverse momentum (pr ~ 2 — 5 GeV/c) in nuclear colli- 
sions, compared with p+p collisions. This observation is well-described by 
models in which baryons and mesons are generated by the recombination 
of quarks drawn from a collectively flowing, thermally equilibrated partonic 
medium. For transverse momenta pr large enough so that hadron masses 
can be neglected, the recombination model predicts hadron spectra to have 



for mesons (M) and baryons (B), respectively. Here w q {wq) denotes the 
phase space density of quarks (antiquarks) in the medium just prior to 
hadronization. These expressions predict meson and baryon yields of ap- 
proximately equal magnitude at the same transverse momentum pr, pre- 
cisely what is observed in Au+Au collisions at RHIC, where the p/n + and 
P/tt~ ratios are found to be near unity at px = 2 — 3 GeV/c. In contrast, 
hadron production by quark or gluon fragmentation predicts baryon-to- 
meson ratios much smaller than unity. 

The recombination model also suggests that the elliptic flow parameters 
of mesons and baryons are related to the elliptic flow parameter of quarks by 
the scaling law (assuming quarks and antiquarks exhibit the same collective 



the form [Ml E3 Ell EZ] : 




flow) [351 [38]: 



«r(Rr)«2^(pr/2); 



4 D V)«3 4 q W3). 
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Indeed, the measured elliptic flow patterns of baryons and mesons show 
remarkable agreement when scaled by the number of valence quarks (see 
Fig. [7|. A combination of hydrodynamics and recombination model con- 
siderations has been used to suggest a slightly modified scaling behavior 
for t>2 when analyzed as a function of transverse kinetic energy E^ — 
(m? +p T ) 1 ^ 2 — m, resulting in two distinct branches, one for mesons and the 
other for baryons (Fig. [TJ left panel). When both V2 and Et are scaled by 
the number of valence quarks (2 for mesons, 3 for baryons) the two branches 
merge into a universal curve for all hadrons (Fig. [7J right panel), indicating 
that the flow pattern is originally developed at the quark levelF] 
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Fig. 7. The transverse kinetic energy (Et = tiit — m) dependences of of various 
mesons and baryons (left panel) collapse into a universal curve when v 2 and Et are 
scaled by the number of valence quarks of the hadron (right panel) . This behavior 
is characteristic for hadron emission by quark recombination from a collectively 
flowing, thermal quark-gluon plasma. 

Jet quenching: QCD jets arise from the hard scattering of incoming 
quarks and gluons and their subsequent fragmentation into directionally 
aligned hadrons. The rates for jet production and other hard-scattering 
processes grow rapidly with increasing collision energy. The production 
rates of such hard probes can be accurately calculated and their interac- 
tions with the medium can be described perturbatively. The access to hard 



1 I am indebted to A. Bialas for the observation that the anisotropy of collective flow 
(elliptic flow) must be manifest at the quark level, if it is established as early in the 
collision (r < 5 fm/c) as hydrodynamics suggests. 
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probes was a primary motivation for constructing RHIC with high center- 
of-mass energy. This strategy has been validated by the discovery of jet 
quenching and its development as a quantitative tomographic probe of the 
quark-gluon plasma. Figure [8] illustrates this discovery, showing the sup- 
pression of 7T° and 7] meson emission in central Au+Au collisions compared 
with expectations from measurements in p+p collisions. The suppression 
of high-pr hadrons stands in distinct contrast to the lack of suppression 
seen in direct photon yields, which are consistent with perturbative QCD 
calculations of their initial production rate. 

In QCD, the energy loss of an highly energetic light parton (quark or 
gluon) is thought to be dominated by gluon radiation induced by the multi- 
ple scattering of the parton on color charges in the medium. The energy loss 
after passage through a homogeneous medium of length L can be expressed 
as AE ~ —a s qL 2 , where the stopping power the medium is governed by 
the parameter [39 

«-'/**«" < 10 > 

Here p denotes the density of color charges in the medium and da/dq 2 is 
the differential scattering cross section for a parton on a color charge, p is 
the inverse color screening length (see eq. ([5])) and Af denotes the mean free 
path of an energetic parton in the medium. The energy loss parameter q 
can also be expressed as a correlation function of gluon fields in the medium 
along the light cone [101 HI] : 

where Cr is the eigenvalue of the SU(3) Casimir operator for the penetrating 
parton. 

The experimental data from RHIC can be used to determine the value 
of q for the quark-gluon plasma produced in nuclear collisions. This de- 
termination requires a realistic modeling of the reaction geometry and its 
time evolution. Such analyses are just now being performed using three- 
dimensional relativistic hydrodynamics |42[ |4"3| |4"4"] . They yield values for 
q in the range of 0.5 — 20 GeV 2 /fm (normalized to the conditions at time 
r = 1 fm/c), considerably larger than the original predictions made on the 
basis of perturbation theory. The large variation of the deduced values of q 
is caused by differences in treating the multiple radiation leading to energy 
loss in the various approaches. 

Additional evidence for the color opacity of the medium is seen in studies 
of the angular correlation of the radiation associated with a high-p-p trigger 
particle |45| 146] . In p+p and d+Au collisions, a hard recoiling hadron 
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Fig. 8. The ratio Raa of measured versus expected yield of various particles 
(7r°,?7,7) in Au+Au collisions at -^/snn = 200 GeV as function of the transverse 
momentum pt- Whereas pions and 77-mesons show the same amount of suppression 
at high pt, where quark fragmentation is the dominant production mechanism, di- 
rect photons are found to be unsuppressed. This indicates that the suppression is 
a final-state effect related to the absorption (energy loss) of energetic partons in 
the medium. 



frequently occurs at 180 degrees in azimuth to the trigger, reflecting the 
back-to-back nature of jets in leading order QCD. In sharp contrast, central 
Au+Au collisions show a strong suppression of such recoils, accompanied 
by an enhancement and broadening of low-pT particle production. Detailed 
analyses indicate that the response of the medium to the passage of an 
energetic parton may be of a characteristic hydrodynamical nature: the 
energy lost by high-energy parton may re-appear as a collective Mach cone. 



4. Toward the "perfect" liquid 

As stated above, a "perfect" liquid is a fluid with the smallest shear 
viscosity allowed by the laws of nature. That quantum mechanics imposes 
a lower limit on the shear viscosity of an entropy carrying fluid can be easily 
seen as follows. In kinetic transport theory, the shear viscosity 77 is governed 
by the rate of momentum transport in the fluid: 



(12) 
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where n is the density, p the average momentum of particles in the medium, 
and <7tr their transport cross section. In quantum mechanics, cross sections 
are bounded by unitarity, which poses a lower bound on the shear viscosity. 
E. g., for s-wave scattering: 



4tt v 3 

°M<~ — *7>~- (13) 
p z 12tt 



A slightly different line of arguments rewrites (12) in terms of the energy 
per particle E/N = e/n and the mean time between scatterings Tf = \{/v, 
and then makes use of the relation S ~ 4iV for the entropy S of a highly 
relativistic system of particles. Finally, the argument recognizes that the 
uncertainty relation sets the bound (E/N)t{ > fi to obtain 

77 Ri — s(E/N)t { > — S S3 — S. (14) 

' 12 v ' ' 12 4vr y ' 

All known materials obey this bound on rj; in fact, until recently, the 



lowest values of r\ for all materials exceeded the bound (14) by at least a 
factor tenj^] Quite generally, the value of r]/s for any material grows at low 
and high temperatures and reaches a minimum somewhere in between. The 
minimum is usually associated with a phase transition; if the transition is 
of first order, a discontinuous jump in r\j s is observed. 

The only system for which the ratio r\js is known to saturate the lower 



bound (14), at least theoretically, is thermal matter governed by the N = 4 
supersymmetric SU(iV c ) Yang-Mills theory in the combined limit of strong 
coupling and large N c [U] . In this limit the theory can be solved exactly by 
mapping it on a weakly coupled superstring theory in five-dimensional Anti- 
deSitter space (more precisely, the space AdSs x S5), the famous AdS/CFT 
duality jH]. The shear viscosity is found to be (setting h = 1) [49]: 

(15) 
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rj 



, 135 C(3) 
+ (8 5 2 iV c )3/ 2 + 



4tt 

There are also more general arguments suggesting that any system that has 



a gravity dual satisfies the bound ( 14 ) [33^ . Quite recently, a real system has 



been found in remarkable parallelism with the developments at RHIC, which 
comes close to the lower bound on rj/s: the dilute, degenerate Fermi gas of 
6 Li atoms at sub-/iK temperatures and artificially induced (via a Feshbach 
resonance) strong coupling among its atoms. The best experimental values 



2 This statement includes helium ( 4 He) below below the lambda point, as long as its 
supcrfluid and normal components flow together in bulk. Although the fraction of 
the entropy carrying normal component tends to zero as T — > 0, its shear viscosity 
increases so rapidly that, overall, n/s — > 00. 
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for this system presently indicate that it comes within a factor of three of 
the rj/s bound. 

What about QCD matter? The shear viscosity of thermal QCD matter 
has been calculated in perturbation theory by evaluating the lowest order 
contributions to parton-parton scattering. At leading logarithmic order and 
for N f = 3 the result is (50]: 

9s 

luvTT a% In a s 

The index "C" is intended to indicate that this is the collisional shear vis- 
cosity. 7/c has also been calculated numerically at leading order a~ 2 beyond 
the leading logarithm |51j . When the result is plotted as a function of T, 
one finds that rj/s should exceed unity for a quark gluon plasma at all tem- 
peratures (see Fig. [9]) [52] . The shear viscosity of a hadronic gas due to 
collisions among pions (and kaons) has also been calculated using, both, ex- 
perimental cross section data|53 and chiral perturbation theory |54| . The 
conclusion is that rj/s may fall slightly below unity at T > 100 MeV but 
rises very rapidly as the temperature decreases (see Fig. [9j. 




10 1 10 s 10 3 10 4 



T(MeV) 



Fig. 9. Predicted shear viscosity-to-entropy density ratio for QCD matter as func- 
tion of temperature. The solid line is the leading-order prediction for a perturbative 
quark-gluon plasma. The dashed line shows the result for a dilute hadronic gas 
using pion cross section data. 
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On the other hand, the RHIC data seem to demand a value of rj/s close 
to the bound (14). When the elliptic flow is calculated in dimensionally 
reduced scenarios - only treating the transverse expansion while maintain- 
ing longitudinal boost invariance - one finds that even a value of r] = s/Att 
seriously affects the agreement with the data obtained in ideal hydrody- 
namics |55| l56l [57] . Obviously, there is a problem. On the one hand, the 
calculation clearly show that the elliptic flow is generated while the mat- 
ter has an energy density significantly exceeding e(T c ) ~ 1 GeV/fm 3 . On 
the other hand, the agreement between data and theory requires rj/s <C 1. 
One concludes that, either the perturbative calculation of the shear viscos- 
ity of a quark-gluon plasma fails by a wide margin, or the shear viscosity 
of the matter produced in nuclear collisions is not dominated by collisions 
among partons. The first alternative has inspired the hypothesis that the 
quark-gluon plasma in the temperature range T c <T < 2T C relevant to the 
nuclear collisions at RHIC is a strongly coupled plasma [58], nicknamed the 
sQGP, which is not composed of partonic quasi-particles and may even not 
be described by quasi-particles at all. 



5. Anomalous viscosity 

The second alternative does not discard the quasi-particle picture alto- 
gether, but to recall that transport processes in electromagnetic plasmas are 
often not dominated by collisions among particles, but by fields permeating 
the plasma. Since the quark-gluon plasma is a plasma, after all, it is not 
inconceivable that similar mechanisms are at work in the matter formed 
at RHIC. Indeed, we have seen earlier (Fig. |4| that plasma instabilities 
generate domains of turbulent color fields whenever the parton momentum 
distribution is anisotropic. This will be the case before thermal equilibra- 
tion is achieved, and it will remain to be true as long as the matter expands 
rapidly in the longitudinal direction. In fact, the momentum anisotropy in 
the later steady state regime is proportional to the shear viscosity. The 
larger the value of r], the more anisotropic is the momentum distribution of 



the expanding quark-gluon plasma (see Fig. 10, right panel). 

Thus, let us assume that coherent domains of color field exist in the 
plasma. The orientation of the field in color space will be random from one 
domain to the next, but the domains will be characterized by an average 
field energy density and a correlation length r m . Let us, for the moment, 
assume that the domain is described by a color-magnetic field B a . The 
momentum deflection of a color charge Q a passing through the domain at 



the speed of light is given by \Ap\ ~ gQ a B a r m , as illustrated in Fig. 10 
(left diagram). Since the domains are assumed to be randomly oriented, 
the momentum direction of a parton will be randomized after traveling a 
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Fig. 10. Left diagram: Thermal partons are deflected when they pass through a 
randomly oriented color field domain of sixe t m . Right diagram: The longitudinal 
expansion in coordinate space implies an oblate momentum distribution of partons 
in the medium. 



distance 

A ( A ) ~ P 2 r ~ t ( 17) 

where again p denotes the average thermal momentum of a parton in the 
plasma. Inserting this expression into the general kinetic theory formula 



(12) for the shear viscosity and substituting s An, p = 3T, we obtain the 



anomalous viscosity [59J: 



9sT 3 



If we use the scaling estimates (B 2 ) ~ g 2 T 4 , r m ~ (gT) 1 mentioned before, 

we obtain tja/s ~ g~ 3 ~ a s 3 ^ 2 , which suggests that the anomalous viscosity 
may dominate over the collisional viscosity in the weak coupling limit. 

In order to obtain a quantitative prediction, we need to calculate the 
anomalous viscosity in the same formal kinetic theory framework as the col- 
lisional viscosity |60| . We start from the perturbed equilibrium distribution 
of partons: 

/(p) = /o(p)[l + /i(p)(l±/o(p))], (19) 

where fo(p) = [exp(n /J p At /T) T 1] _1 is the Bose (gluons) or Fermi (quarks) 
distribution. The perturbation fi(p) for shear strain is parametrized by the 
shear viscosity n: 

fi(p) = (ViU k + V kUi - p lk V • u) (20) 
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The phase space distribution of partons, f(p), is governed by the non-abelian 
Vlasov-Boltzmann equation 

[d t + vV r + F- V p ] /(r, p, t) = C[f] , (21) 

where v = p/E p is the parton velocity, F = gQ a {E a + v x B a ) the color 
Lorentz force, and C[f] denotes the collision term. If we assume that the 
color fields are randomly oriented and distributed, the ensemble averaged 
phase space distribution f{p) can be shown to satisfy an equation of the 
Fokker-Planck-Boltzmann type: 

[d t + v ■ V r - V p ■ D(p) ■ V p ] fir, p, t) = C[f] , (22) 

where the momentum diffusion tensor Dik is given by 

Afc(p)=/* dt'(Fir',t')F k ir,t)) with r' = r + vtf - t). (23) 



For isotropically distributed color field domains with a correlation length r m 
the diffusion tensor takes the form Di^ip) = (F 2 )r m . Following the standard 
Chapman-Enskog approach for the calculation of transport coefficients, we 
obtain the following result for the shear viscosity: 

v, = 0(1) (4 -i)^ + 0(10 ' s ^a + W (24) 

This expression exhibits several remarkable features. First, the numer- 
ator of the anomalous contribution to the shear viscosity, {F 2 )r m / (iV^ — 1), 
is precisely the contribution of the color field domains to the energy loss 
parameter q. This is easily seen from the formulation of q as a gluon cor- 
relation function along the light cone, eq. ( |11[ ). Secondly, a more detailed 
scaling analysis reveals that the quantities (F 2 ) and r m are functions of the 



momentum space anisotropy of the plasma and thus, by virtue of (20), of 
r\ itself. Taking this implicit dependence into account, one finds that the 
anomalous viscosity scales as 

rj A /s ~ ig 2 \Vu\)-^, (25) 

where | Vit| denotes the magnitude of the shear strain. The explicit depen- 
dence of t]a on the shear strain gives meaning to the term "anomalous" vis- 
cosity. Since the normal, collisional viscosity scales as rjc ~ (s ,4 ln g^ 1 )^ 1 , 
but is independent of |Vit|, the collisional viscosity always dominates for 
very small shear strains, whereas the anomalous viscosity dominates for 
large shear strains and weak coupling! More quantitative statements re- 



quire the numerical evaluation of the correlator (23) 
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6. Exploring the "perfect" liquid 

It turns out that the relationship between q and ry noticed in the case 
of the anomalous viscosity can be generalized to any situation where the 
thermal excitations of the medium have the same quantum numbers and 
interactions as highly energetic ("hard") excitations. It not only applies to 
QCD, but to any unbroken gauge theory, in which transport cross sections 
are dominated by small angle scattering. Under these conditions, one can 
derive the following relation |61j : 



1 
s 



1.25 



(26) 



This correspondence fails to hold when the medium is so strongly coupled 
that its thermal excitations cannot be described as quasi-particles, as in 
strongly coupled super- Yang-Mills (SYM) theory, or when the thermal ex- 
citations are described by quasi-particles of a different kind, as in QCD at 



temperatures below T c . This is illustrated in Fig. 11 where the left-hand 



side and the right-hand side of eq. (26) are plotted separately versus the 
temperature (for QCD) or the 't Hooft coupling A = g 2 N c (for the SYM 
theory). In both cases, the two sides of the equation diverge when the 
medium is not described by a partonic quasi-particle perturbation theory. 



In both instances, the right-hand side of (26) is a better measure of the 



fundamental coupling strength when the relationship fails to hold. 





QCD 




/(in ry 4 







1 



Mr/r e ) 




]n(l/A) 



Fig. 11. The left-hand side (rj/s, red lines) and right-hand side (1.25 T 3 /q, green 



lines) of eq. (26 1 plotted versus the effective coupling strength. Left panel: QCD; 
right panel: N = 4 super- Yang-Mills theory. For QCD, the abscissa shows T/T c 
on a logarithmic scale; for the SYM theory, the abscissa shows the inverse 't Hooft 
coupling, also on a logarithmic scale. 
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Reliable independent determinations of rj/s and q/T 3 from the RHIC 
data would thus permit a model independent determination of the strongly 
coupled nature of the matter produced in the collisions. If the relation 
(26) holds, the quark-gluon plasma at RHIC is composed of partonic quasi- 
particles; if it fails, it is a strongly coupled quark-gluon plasma (sQGP). 
Presently, the analysis of the RHIC data in terms of rj/s and q/T 3 is not 
yet precise enough to make this determination possible. Using the value go ~ 
1—2 GeV 2 /fm extracted in the higher-twist analysis of jet quenching[62] and 
the initial temperature To ~ 335 MeV deduced from the measured entropy 
of the final state[63], one finds 1.25 T 3 /q « 0.12 — 0.24, not inconsistent with 
the constraints on rj/s obtained from comparisons of viscous hydrodynamics 
calculations with the elliptic flow data. 

If the effect of the medium on a hard parton can tell us something about 
the properties of the matter formed at RHIC, then the fate of the energy lost 
by the parton must be even more informative. This issue has two aspects: 
the energy lost by the triggered parton and the energy lost by its recoil 
partner. The latter is far larger than the former because, on average, the 
triggered parton originates close to the near-side surface of the fireball and 
thus traverses only a relatively small amount of matter, while the recoil 
partner, again on average, propagates through the bulk of the fireball. In 
fact, the recoil parton loses so much of its energy that its fragments become 
virtually indistinguishable from the matter itself. However, the energy does 
not disappear, and its kinematic distribution can be identified by careful 
background subtraction. The surprising result of this analysis was that 
the excess energy does not appear at 180° in azimuth with respect to the 
trigger, but at an angle of about 110° (see Fig. 12 A three-body coincidence 
analysis indicates that the energy appears on both sides of the recoil, ruling 
out a simple deflection mechanism. 

The explanation most compatible with the data, but also the most spec- 
ulative and spectacular one, is that the energy lost by the recoil parton 
travels through the medium in the form of a sonic Mach cone [64 . If this is 
true, the angular peak position can be used to determine the sound velocity 
of the quark-gluon plasma! It would also show that sound propagation in 
the medium is only weakly damped and thus confirm that the shear vis- 
cosity is small. The formation of a Mach cone trailing a fast color charge 
has been demonstrated by explicit calculations in the strongly coupled SYM 
theory [65]. In QCD, the precise form of the coupling between a fast color 
charge and the sound mode has not yet been derived. What is clear from 
phenomenological studies is that the coupling must be highly (about 70%) 
effective in terms of energy transfer [66 . 

The emission pattern of secondaries on the trigger side is also modified 
in Au+Au collisions compared with p+p collisions. In addition to a jet-like 
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Fig. 12. Angular distribution of secondaries around the beam axis with respect 
to the direction of the trigger hadron. The red part shows the jet accompanying 
the triggered hadron, the blue part shows the fragments of the recoil parton. Left 
panel: p+p collisions exhibit the usual 180° di-jet event structure. Right panel: In 
Au+Au collisions the peak of the recoil secondaries is shifted to 110° in azimuth, 
suggestive of a Mach cone around the recoil direction. 



pattern of secondaries at small angles with respect to the trigger hadron with 
nearly the same characteristics as in p+p collisions, one finds also a wide 
ridge near 0° in azimuth but distributed over a wide pseudorapidity range 
(A?7 > 3) [67]. The chemical composition and spectral distribution of the 
hadrons making up the ridge is very similar to the bulk matter. One possible 
explanation is that the ridge is formed by the energy lost by the parton 
that produced the trigger hadron. Since the jet cone itself has the same 
shape in p+p and Au+Au, it must be formed in vacuum, after the parton 
escapes from the fireball. But why does the energy radiated during its 
passage through the matter spread out longitudinally, but not in azimuth? 
This could actually be the result of the deflection of the radiated gluons by 
turbulent color fields present in the quark-gluon plasma. If these are mostly 
color-magnetic fields oriented in the transverse direction to the beam axis, 
as it is suggested by a perturbative analysis of the color field instabilities, 
they would result in a preferentially longitudinal deflection of the radiated 



gluons (see Fig. 13, left panel). The diffusion in pseudorapidity rj, but 



not in azimuthal angle (j) is, indeed, found when one solves the diffusive 
transport equation ( 22 ) for radiated gluons in the presence of transversely 
oriented color-magnetic fields (see Fig. 13 right panel) [68]. Similar effects 
have recently also been found in self-consistent microscopic simulations of 
partons propagating in spontaneously generated color fields [69 . 
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Fig. 13. Left panel: Schematic representation of the longitudinal deflection of radi- 
ated gluons propagating through transversely oriented color-magnetic fields. Right 
panel: Longitudinal broadening of the distribution of radiated gluons obtained by 
solution of the diffusive Boltzmann equation ( 22 1 . The final pseudorapidity distri- 
bution is shown by the solid red curve. 



7. Conclusion 



The incredible wealth and quality of the RHIC data have challenged the- 
oretical ideas of the physical properties of the quark-gluon plasma on many 
fronts. The most remarkable insight gleaned from the data is that the mat- 
ter produced in nuclear collisions at RHIC is an almost " perfect" , i. e. nearly 
inviscid liquid, whose shear viscosity approaches the lower bound dictated 
by unitarity. The strongly coupled super- Yang- Mills theory provides one 
model for such a liquid; an alternative model is the turbulent plasma per- 
meated by color fields generated by expansion driven plasma instabilities. 
If both, shear viscosity rj and energy loss parameter q of the matter can 
be determined reliably from an analysis of the data, it may be possible to 
distinguish between the two scenarios. The found relation, connecting a 
large color opacity of the quark-gluon plasma with a low dissipation for col- 
lective flow, is clearly confirmed by the RHIC data. Increasingly differential 
studies of the phenomena accompanying jet formation in nuclear collisions 
are promising further surprises, which will keep theorists occupied, if not 
puzzled, for years to come. 
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